The interleukin-6/interleukin-23/T helper 17-axis as a driver of neuro-immune toxicity in the major neurocognitive psychosis or deficit schizophrenia: A precision nomothetic psychiatry analysis

Background Schizophrenia and especially deficit schizophrenia (DSCZ) are characterized by increased activity of neuroimmunotoxic pathways and a generalized cognitive decline (G-CoDe). There is no data on whether the interleukin (IL)-6/IL-23/T helper 17 (IL-6/IL-23/Th17)-axis is more associated with DSCZ than with non-deficit schizophrenia (NDSCZ) and whether changes in this axis are associated with the G-CoDe and the phenome (a factor extracted from all symptom domains) of schizophrenia. Methods This study included 45 DSCZ and 45 NDSCZ patients and 40 controls and delineated whether the IL-6/IL-23/Th17 axis, trace elements (copper, zinc) and ions (magnesium, calcium) are associated with DSCZ, the G-CoDe and the schizophrenia phenome. Results Increased plasma IL-23 and IL-6 levels were associated with Th17 upregulation, assessed as a latent vector (LV) extracted from IL-17, IL-21, IL-22, and TNF-α. The IL-6/IL-23/Th17-axis score, as assessed by an LV extracted from IL-23, IL-6, and the Th17 LV, was significantly higher in DSCZ than in NDSCZ and controls. We discovered that 70.7% of the variance in the phenome was explained by the IL-6/IL-23/Th17-axis (positively) and the G-CoDe and IL-10 (both inversely); and that 54.6% of the variance in the G-CoDe was explained by the IL-6/IL-23/Th17 scores (inversely) and magnesium, copper, calcium, and zinc (all positively). Conclusion The pathogenic IL-6/IL-23/Th17-axis contributes to the generalized neurocognitive deficit and the phenome of schizophrenia, especially that of DSCZ, due to its key role in peripheral inflammation and neuroinflammation and its consequent immunotoxic effects on neuronal circuits. These clinical impairments are more prominent in subjects with lowered IL-10, magnesium, calcium, and zinc.

There are also some reports that trace elements such as copper and zinc and ions such as calcium and magnesium, as well as serum albumin, may play a role in schizophrenia. Thus, serum copper is elevated in schizophrenia patients [41], while lower zinc is shown in a recent meta-analysis [42]. Patients with schizophrenia often have lowered serum and intracellular magnesium levels, while magnesium treatment may improve psychosis [43]. There is also some evidence that aberrations in calcium signaling are associated with neurocognitive impairments [44] and that lowered albumin may contribute to the severity of the phenome of schizophrenia [45].
Nevertheless, there is no data on whether the IL-6, IL-23, and Th17 (IL-17, TNF-α, IL-21 and IL-22) axis is more associated with MNP than with SNP and whether changes in this axis are associated with the G-CoDe and the phenome of schizophrenia. Hence, this study was conducted to delineate a) whether the IL-6/IL-23/Th17 axis is activated in MNP and is associated with impairments in the G-CoDe and the phenome; b) the associations between IL-6 and IL-23 and Th17 effector cytokines; and c) the role of IL-10, copper, zinc, calcium and magnesium in the neuro-immune pathophysiology of MNP.

Participants
Ninety patients with schizophrenia and forty healthy controls were included in this study. Patients were recruited between February and June 2021 at "The Psychiatry Unit," Al-Hakeem General Hospital, Najaf Governorate, Iraq. All patients were stable for at least 12 weeks. Family and friends of personnel or friends of patients served as controls. Patients and controls were excluded if they had ever taken immunosuppressive treatments or glucocorticoids or if they had been diagnosed with a neurodegenerative or neuroinflammatory illness such as Alzheimer's disease, Parkinson's disease, multiple sclerosis, or stroke. Additionally, individuals with (auto)immune diseases such as inflammatory bowel disease, rheumatoid arthritis, COPD, psoriasis, or diabetes mellitus were excluded. None of the controls had a present or lifetime DSM-IV-TR axis I diagnosis or a family history of schizophrenia or psychosis. Patients having axis-1 DSM-IV-TR disorders such as bipolar disorder, major depression, schizoaffective disorder, obsessive-compulsive disorder, psycho-organic illnesses, or substance use disorders were excluded. CRP levels were < 6 mg/L in all subjects, indicating that there was no overt inflammation.
The study followed Iraqi and international privacy and ethics laws. Before participating in this study, all participants and first-degree relatives of participants with schizophrenia gave written informed consent (legal representatives are mother, father, brother, spouse, or son). The study was approved by the ethics committee (IRB) of the College of Science, University of Kufa, Iraq (82/2020), which follows the Declaration of Helsinki's International Guideline for Human Research Protection.

Clinical assessments
To collect patient and control data, a senior psychiatrist with expertise in schizophrenia conducted semi-structured interviews. Using DSM-IV-TR diagnostic criteria, the Mini-International Neuropsychiatric Interview (M.I.N.I.) was used to diagnose schizophrenia. The same psychiatrist evaluated the Schedule for the Deficit Syndrome (SDS) [46] for the diagnosis of primary deficit schizophrenia, as well as the Positive and Negative Syndrome Scale (PANSS), the Scale for the Assessments of Negative Symptoms (SANS) [47], the Brief Psychiatric Rating Scale (BPRS) [48], and the Hamilton Depression (HAM-D) and Anxiety (HAM-A) rating scales [49,50]. We calculated z-unit weighted composite scores reflecting psychosis, hostility, excitation, mannerism, PMR and FTD symptom domains as previously reported using BPRS, HAM-D, HAM-A, and PANSS items and the total SANS score was used to assess negative symptoms [7,8,23,51]. Based on the 6 computed z composite scores reflecting PHEM, negative, PMR, and FTD scores, we then extracted the first factor from these data, named the "phenome". On the same day, a research psychologist performed neuropsychological probes using the Brief Assessment of Cognition in Schizophrenia (BACS) [52] while remaining blind to the clinical diagnosis. The latter battery consists of the Digit Sequencing Task (which measures working memory), the List Learning test (which tests verbal episodic memory), Controlled Word Association (which assesses letter fluency), Category Instances (which measures semantic fluency), the Tower of London (to probe executive functions), Symbol Coding (which probes attention), and the token motor task. The total score was computed. Based on our previous findings [5] that one latent vector underpins all cognitive tests, we here extract the first factor from the 7 BACS test, and named this factor the "G-CoDe" (generalized cognitive decline). Tobacco use disorder (TUD) was diagnosed in accordance with DSM-IV-TR criteria. The following formula was used to compute the body mass index (BMI): body weight (kg) / length (m 2 ).

Assays
Each subject's fasting venous blood was obtained in the early morning hours. After 15 minutes at room temperature, the blood was allowed to coagulate for 10 minutes before being centrifuged at 3000 rpm for 10 minutes. Separated serum was then transferred to Eppendorf tubes and stored at -80˚C until analysis. The concentrations of CRP in serum were measured using a kit supplied by Spinreact 1 (Barcelona, Spain). The test is based on the notion of latex agglutination. Serum albumin, calcium, and magnesium concentrations were determined using spectrophotometric kits provided by Biolabo 1 (Maizy, France). Copper was determined spectrophotometrically by kits supplied by LTA Co., (Milano, Italy). Zinc concentrations were determined using a kit provided by Cenrionic GmbH (Wartenberg, Germany). We measured serum IL-6, IL-10, G-CSF, and IL-1β using commercial ELISA sandwich kits provided by Sunlong Biotech Co., Ltd. (Zhejiang, China), and Melsin Medical Co. (Jilin, China) provided the  other ELISA kits (IL-17, IL-21, IL-22, IL-23, and TNF-α). We diluted samples as required in the presence of very concentrated biomarkers. The within-assay coefficient of variation for all assays was <10.0% (precision within-assay). The sensitivities of the kits were 1.0 pg/ml for IL-10, IL-17, IL-21, IL-22, IL-23, TNF-α, G-CSF, and 0.1 pg/ml for IL-1β, and IL-6.

Statistical analysis
To compare scale variables between groups, one-way analysis of variance was employed, while analysis of contingency tables (two tests) was used to analyze category variables. We used multivariate general linear model (GLM) analysis to examine the association between the 14 biomarkers and the diagnosis while entering age, gender, education, BMI, TUD and drug use as covariates. The between-subject effects tests were employed to investigate the impact of significant explanatory factors on each of the observed biomarkers. As a result, we calculated the estimated marginal means (SE) values produced by the model (adjusted for the significant confounders). These multiple associations were subjected to p-correction for false discovery rate (FDR) [53], and the protected Least Significant Difference (LSD) tests were used to assess multiple pair-wise differences. Multiple regression analysis (automated approach with p-toentry of 0.05 and p-to-remove of 0.06 while assessing the change in R 2 ) was used to determine the important biomarkers that predict the phenome in schizophrenia and all participants combined. Multivariate normality (Cook's distance and leverage), multicollinearity (using tolerance and VIF), and homoscedasticity (using White and modified Breusch-Pagan tests for homoscedasticity) were checked. These regression analyses' results were always bootstrapped using 5.000 bootstrap samples, and the latter results are shown if the results were not concordant. We used two-step cluster analysis to derive meaningful clusters from the data set. Feature reduction was performed using principal component (PC) analysis. IBM SPSS Windows version 25, 2017 was used to conduct the statistical analysis.
Partial Least Squares (PLS)-SEM pathway analysis is a statistical approach for predicting complex cause-effect linkages utilizing both discrete indicators (variables) and latent variables (factors based on a collection of highly connected indicators) [54]. Without imposing distributional assumptions on the data, PLS allows the estimation of complex multi-step mediation models with several latent constructs, single indicator variables, and structural pathways (associations between indicators or latent vectors). This strategy has recently been used to develop novel nomothetic models of affective disorders and schizophrenia by combining the many components of a disease into a causal-effect, mediation model [9,22]. A causative framework is developed utilizing causome, protectome, adverse outcome pathways, and phenome indicators (either single indicators or latent vectors), which is then analyzed and cross-validated using PLS pathway analysis on bootstrapped data (e.g., 5.000 samples). As a result, pathway coefficients (with exact p values), specific indirect (mediated) and total effects are estimated to determine the influence of direct and mediated pathways. According to the power analysis, the predicted sample size for multiple regression analysis (which applies to PLS) should be at least 70 to obtain a power of 0.8, an effect size of 0.2, and an alpha of 0.05 with 5 preset variables. Complete path analysis with 5.000 bootstrap samples was conducted only when the outer and inner models matched the following quality criteria: a) model quality as measured by the SRMR index is less than 0.08; b) outer model loadings on the latent vectors are greater than 0.666 at p < 0.001; and c) the latent vectors exhibit accurate construct validity as measured by average variance extracted > 0.5, Cronbach's alpha > 0.7, rho A > 0.7, and composite reliability > 0.8. Compositional invariance was investigated using Predicted-Oriented Segmentation analysis, Multi-Group Analysis, and Measurement Invariance Assessment. The predictive performance of the model was evaluated using PLSpredict with a 10-fold cross-validation. Table 1 shows the sociodemographic and clinical data of MNP versus SNP and controls. There were no significant differences in age, sex, BMI, marital status, TUD, residency, and employment status between MNP and SNP and healthy controls. Both patient groups were somewhat less educated than the control group, while there were no significant differences in education level between both patient groups. MNP patients have a significantly higher (p<0.001) positive family history of schizophrenia than the SNP group. Since some of the patients were treated with olanzapine (n = 62), clozapine (n = 12), fluphenazine (n = 27), haloperidol (n = 7), quetiapine (n = 5), risperdal (n = 16), benzodiazepines (n = 19), and mood stabilizers (n = 39) we have examined the impact of these drugs on the clinical scores. After adjusting for age, sex, BMI, smoking, education, and the drug state of the subjects, the major symptom domain scores (psychosis, hostility, excitement, mannerism, SANS, PMR, and FTD) were significantly different between the three study samples and increased from controls ! SNP ! MNP. These differences remained significant after FDR p-correction at p<0.001.

Differences in biomarkers between the study groups
Multivariate GLM analysis ( Table 2) was used to assess the relationships between biomarkers and diagnosis after controlling for age, BMI, sex, smoking and drug state. There were highly

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia significant differences in the biomarkers across the groups (which remained significant after FDR p-correction at p<0.012), but no significant impact of the covariates.Tests for betweensubject effects in Table 3, which shows the estimated marginal means, revealed that IL-23, IL-6, IL-17, and TNF-α were significantly higher in MNP than in SNP and controls and that IL-17 and TNF-α were higher in SNP than in controls. IL-1β, IL-22, G-CSF, IL-21, and IL-10 were significantly higher in schizophrenia than in controls. There were no differences in albumin between the study groups, and zinc levels were lowest in both schizophrenia groups. Magnesium was significantly lower in MNP as compared with the other two groups, whereas calcium was lower in MNP than in SNP. Copper levels were significantly higher in MNP than in SNP. We found no significant effects of the drug state on the single biomarkers (tested with multivariate GLM analysis and tests for between-subjects effects) even without p correction for FDR. In this model, IL-23 and IL-6 were considered to be correlated with a latent vector (LV) extracted from IL-17, IL-21, IL-22, and TNF-α (labeled Th17 LV) while the latter predicted G-CSF and IL-10. After feature selection, multi-group analysis, PLS predict analysis, and The results of these GLM analyses were adjusted for use of olanzapine, clozapine, modecate, kemadrine, risperidone, benzodiazepines, mood stabilizers (not shown as all non-significant).

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia prediction-oriented segmentation, Fig 1 shows the results of the PLS path analysis conducted on 5.000 bootstrap samples. The final PLS model did not contain non-significant pathways and variables. With SRMR = 0.039, this model's overall fit was more than sufficient. Additionally, the reliability of the factor constructs was satisfactory for all three factors, with Cronbach alpha > 0.754, rho A > 0.785, composite reliability > 0.845, and AVE > 0.58. At p <0.0001, all loadings on the outer models were more than 0.620. We discovered that 70.7% of the variance in the phenome could be explained by regression on the cognitome, Th17 LV (both positively) and IL-10 (inversely), and that 68.9% of the variance in the cognitome could be explained by education (positively) and the combination of IL-23 and Th17 LV. IL-23 and IL-6 explained 17.7% of the variance in the Th17 LV. The cognitome (0.558), the Th17 LV (0.093), and the phenome (0.624) all had acceptable construct cross-validated redundancies (as determined by PLS blindfolding). Confirmatory Tetrad Analaysis (CTA) in PLS demonstrated that the latent vectors were not conceptualized incorrectly as reflective models. PLS prediction revealed that the indicators of the endogenous constructs had positive Q 2 predict values, indicating that they outperformed the naïve benchmark. There were significant specific indirect effects of IL-23 on the phenome through the Th17 LV (t = 2.82, p = 0.002), the cognitome (t = 1.68, p = 0.047), the path from Th17 LV to the cognitome (t = 3.64, p<0.001), and the path from Th17 LV to IL-10 (t = -1.84, p = 0.033). Although the direct pathway from IL-6 to the phenome was insignificant, the overall impact of IL-6 on the phenome was significant (t = 1.74, p = 0.041). Nonetheless, IL-23 (t = 6.76, p<0.001) and the Th17 LV (t = 5.10, p<0.001) had a much greater influence on the phenome.

Construction of pathway endophenotype and endophenotype classes
Based on the PLS findings that IL-6 and IL-23 are both associated with the Th17 LV we constructed an LV extracted from IL-6 and IL-23 (AVE = 58.7%, composite reliability = 0.738,

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia loadings: 0.823 and 0.705, respectively) and then extracted an LV from this construct and the Th17 LV, yielding a new IL-6/IL-23/Th17 LV (AVE = 70.8%, loadings: 0.841 and 0.841, respectively, composite reliability = 0.844). Tables 2 and 3 show that this IL-6/IL-23/Th17 score significantly differed between the three study groups and increased from controls ! MNP ! SNP. In addition, we were able to extract a latent vector from the IL-6/IL-23/Th17 (loadings: 0.880), G-CoDe (-0.918) and phenome (0.918) scores (AVE = 81.52%, composite reliability = 0.912). Two-step cluster analysis with the IL-6/IL-23/Th17, G-CoDe, and phenome entered as continuous variables and the diagnosis of schizophrenia as a categorical variable showed three clusters of patients with a silhouette measure of cohesion and separation of 0.66. The formed schizophrenia clusters exactly matched with the MNP (deficit) versus SNP (nondeficit) dichotomy. We found no significant effects of these drugs on the IL-6/IL-23/Th17 score (tested with univariate GLM analysis).

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia

Prediction of symptom domains and G-CoDe by the IL-6/IL-23/Th17 axis score
The first multiple regression analysis in Table 4 entered the phenome LV as the dependent variable and the IL-6/IL-23/Th17 LV and G-CoDe scores as the explanatory variables while allowing for the effects of sex, age, smoking, education, the trace elements, and ions. Table 4, regression #1 demonstrates that 68.8% of the variance in the phenome was explained by the regression on G-CoDe and IL-6/IL-23/Th17 scores . Figs 2 and 3 show the partial regression plots of the phenome on both predictors (after considering the effects of the other covariate).
In the restricted study sample of schizophrenia patients (regression #2), the same variables explained 28.5% of the variance in the phenome. Table 4, regression #3 demonstrates that 54.6% of the variance in the G-CoDE was explained by the regression on IL-6/IL-23/Th17 (inversely) and magnesium, copper, and zinc (all positively). In the restricted study sample of schizophrenia patients (regression #4), we found that the IL-6/IL-23/Th17 score (inversely) and magnesium and calcium (positively) explained 28.5% of the variance in the G-CoDe.

Increased cytokine levels in MNP
The first major finding of this study is that IL-23, IL-6, IL-17, and TNF-α were considerably greater in MNP than in SNP, while IL-1β, IL-22, G-CSF, IL-21, IL-17, IL-10, and TNF-α were significantly higher in schizophrenia than in controls. Furthermore, elevated IL-23 and IL-6 were identified in MNP but not SNP as compared to controls, suggesting that these cytokines are characteristics of MNP rather than schizophrenia in general. Previous research has linked

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia MNP to increased TNF-α, IL-6, and IL-1β production, as well as CCL11, CCL2, and the soluble IL-1 receptor antagonist [17, 55, 56], and increased TNF-α levels have been identified as a critical factor in the pathophysiology of MNP [57]. Increments in cytokines such as IL-1β, TNF-α, IL-6, and IL-10 have frequently been described in schizophrenia [58]. Blood levels of IL-23 were found to be high in people with treated schizophrenia, FEP, and repeated episodes of schizophrenia [28, 59-61]. As reviewed by Roomruangwon et al. [58] and reported in some [27, 28, 30, 60, 62], but not all [63], studies, schizophrenia is accompanied by increased IL-17 levels. Increased IL-22 and IL-6, but not IL-17 or IL-1β, were observed in schizophrenia as compared with healthy subjects [64]. Interestingly, antipsychotic treatment for 3 months was found to significantly reduce plasma levels of Th17-pathway-related cytokines like IL-6 and IL-17A [65]. Contradictorily, a four-week treatment with antipsychotic agents did not affect circulating levels of IL-23 or IL-17 [66,67]. Most likely, the differences mentioned above come from comparing different research groups, like MNP vs. SNP/controls or schizophrenia vs. controls.

IL-23 and IL-6 as drivers of Th17 pathway activation
The second major finding is that our PLS analysis revealed that increasing IL-23 and IL-6 levels are linked to Th17 activation, as measured by extracting a latent vector from IL-17, IL-21, IL-

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia 22, and TNF-α. As described in the introduction, IL-23 and IL-6 are key factors in the growth and survival of the pathogenic Th17 phenotype [34-37, 67,68]. Moreover, IL-6 and IL-23 signaling is required to activate STAT3 and IL-17 production, while IL-21 activates STAT3 [69,70]. IL-23 is a key cytokine involved in a variety of inflammatory and autoimmune disorders by producing a highly pathogenic T cell population [71]. Th17 cells are the main producers of IL-22 [72] and IL-17 coupled with IL-22 marks a particularly pathogenic population of aggressive autoreactive and proinflammatory T cells [73]. Nevertheless, IL-22 also has protective properties [74], including acting as an antimicrobial peptide and collaborating with IL-17 in regulating the immune response [72,73]. IL-21 is another proinflammatory cytokine that mediates antibody class switching and production, and this cytokine has pathogenic effects in (auto)immune disorders, including Sjogren's syndrome, systemic lupus erythematosus, and psoriasis [75][76][77].
As a result, rather than calling this axis the TNF/IL-23/IL-17 axis [78] or the IL-23/IL-17/ G-CSF axis [66], it is more appropriate to refer to it as the pathogenic IL-6/IL-23/IL-21/IL-17/ IL-22/TNF-axis, or simply the IL-6/IL-23/Th17-axis. In this way, we were able to find a common core between IL-23, IL-6, and the Th17 cytokines, which confirmed the concept.
It is worth noting that in schizophrenia, activation of the pathogenic IL-6/IL-23/Th17 axis is linked to other important pathways. First, brain-derived neurotrophic factor (BDNF) is adversely related to IL-23 and illness severity ratings in schizophrenia [79], and a first proteinprotein interaction (PPI) network created using FEP/FES genes shows that BDNF is part of the same immune PPI network as STAT3, IL-6, TNF-α, IL-1β, and IL-10 [80]. Second, complement component factors that have a role in schizophrenia, such as C1q, C3, and C4 [21, 81,82], may affect the production of IL-23 and IL-17 family members [83]. Third, blood-brain

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia barrier (BBB) endothelial cells express IL-17 and IL-22 receptors, and binding to IL-17 and IL-22 causes BBB permeabilization [84,85], thereby aggravating the effects of lipopolysaccharides (LPS), other cytokines and tryptophan catabolites on BBB breakdown [86]. Fourth, the IL-6/ IL-23/Th17 axis plays an important role in the immune-inflammatory responses directed to gut commensal microbiota (including Gram-negative bacteria). Thus, commensal bacteria may prime IL-23 expression and production by dendritic cells [87], and induce IL-22 production in the gut [88]. IL-23 increases gut permeability by a breakdown of the tight junction barrier [89] and orchestrate the immune-inflammatory response in the gut [90]. IL-22 generally has more regenerative and protective functions [91][92][93], including the regulation of antimicrobial peptide production and the composition of gut commensal bacteria [88]. IL-17 plays a key role in response to infectious agents by mounting an innate immune response in epithelial cells, which may lead to systemic inflammatory and autoimmune responses [91]. IL-21 may trigger gut inflammation and plays a role in inflammatory bowel disease, while claudin-5 is a downstream gene of this cytokine [94]. Breakdown of the tight junctions and adherents barriers (leaky gut) with increased bacterial translocation is involved in schizophrenia and, especially, in MNP [86,95].

Effects of the pathogenic IL-6/IL-23/Th17 axis on G-CoDe and the phenome of schizophrenia
The third major finding of this study is that the latent vector score extracted from the IL-6/IL-23/Th17 data was considerably greater in MNP than SNP, with a difference of more than one

PLOS ONE
IL-6/IL-23/Th17-axis as a driver of neuroimmune toxicity in the MNP or deficit schizophrenia standard deviation, and predicted cognitive deficits and the severity of the phenotype of schizophrenia. Moreover, we showed that a reliable latent vector could be extracted from the IL-6/IL-23/Th17, G-CoDe and phenome scores, thereby shaping a new pathway phenotype reflecting the link among peripheral adverse outcome pathways and the cognitome and phenome of schizophrenia. Furthermore, clustering analysis revealed two qualitatively distinct and well-defined patient groups using the same three scores, and these groups perfectly matched the deficit vs. non-deficient dichotomy. These findings support the view that MNP and SNP are qualitatively distinct classes [96] and that the pathogenic IL-6/IL-23/Th17 axis is a major player in this regard.
We previously addressed IL-6 and TNF-'s many neurotoxic actions, which, when coupled with other neurotoxic substances like LPS and lipid peroxidation, may produce neurotoxicity and, hence, explain the cognitive deficits and phenome of schizophrenia [57,97]. As discussed previously, these compounds may induce a plethora of neurotoxic effects, especially on processes like neurogenesis, neuroplasticity, cerebral cortex radial glia guided migration, synapse assembly, axogenesis, axonal spreading and branching, synapse assembly and structure, preand post-synaptic neuronal connectivity, regulation of excitatory synaptic functions, and postsynaptic protein assembly [21]. By inference, the strong effects of the IL-6/IL-23/Th17 axis on both the cognitive detrioration and symptomatome of schizophrenia indicate that IL-23, Th17, IL-22 and IL-21 may further contribute to the neurotoxic effects.
A variety of cells in the central nervous system (CNS) produce IL-17 under physiological conditions, while IL-17R is expressed by immune cells and glial cell populations [98]. IL-23 is

PLOS ONE
produced by CNS cells and is a principal mediator of microglial activation, neuroinflammation, and tissue damage [99]. There are two ways in which Th17 cells are involved in neuroinflammation, either directly via the production of cytokines (IL-17, IL-21; IL-22; IFN-γ; and G-CSF) or indirectly via the stimulation of neutrophil infiltration and microglial cells, which secrete cytokines, and by attracting CD8+ and Th1 cells to the CNS [25]. IL-17 has been hypothesized to adversely influence adult hippocampus neurogenesis even in physiological conditions [100]. Moreover, IL-23-activated Th17 cells in the brain are essential for sustaining persistent neuroinflammation during infection and autoimmune responses [34]. Furthermore, following stroke, IL-21 and its receptor are expressed in the CNS via brain infiltrating CD4+ T cells, and CNS IL-21 strongly contributes to CNS tissue damage [101]. IL-22 is expressed in the brain in physiological conditions, and is upregulated during neuroinflammatory conditions [102,103]. While this cytokine may have protective effects [103,104] and induces an acute phase response [105], it may also display more detrimental inflammatory effects [88]. Most importantly, during inflammatory responses, IL-22 is upregulated in the brain and may increase the production of TNF-α, IL-6, and prostaglandins and induce STAT3, MAP-kinase, and JAK-STAT pathways [102], which all play a role in schizophrenia [102,106].
It is also important to note that after considering the detrimental effects of the IL-6/IL-23/ Th17 axis on the phenome of schizophrenia, IL-10 showed an inverse association with the phenome. Since IL-10 levels were higher in patients than in controls and since IL-10 is a negative immunoregulatory cytokine [16], our results indicate that a relative shortage in IL-10 CIRS functions may aggravate the effects of the pathogenic IL-6/IL-23/Th17 axis.

Calcium, magnesium and zinc in schizophrenia
This study shows that schizophrenia is accompanied by lower zinc as compared with controls, and that MNP is characterized by lower magnesium and calcium but higher copper as compared with SNP. As described in the Introduction, some [42, [107][108][109][110][111][112][113] but not all [114,115] authors reported changes in these trace elements and ions in plasma, serum or CSF in schizophrenia patients. Although some studies reported inverse associations between calcium and magnesium levels and the severity of schizophrenia, others could not find any association [107,116]. Still, we found that magnesium, zinc, copper, and calcium were positively linked to the G-CoDe. This means that these trace elements and ions protect cognitive functions from the neurotoxic effects of the IL-6/IL-23/Th17 axis, and that the lower levels in schizophrenia and MNP contribute to the overall decline in cognitive function.
It is important to note that inflammation or the acute phase response is accompanied by lowered levels of zinc, magnesium, and calcium (reviews in: [117,118]. Moreover, these three elements and copper all modulate NMDA receptor (NMDAR) functions: a) zinc and magnesium are anatagonists of the glutamatergic NMDAR [119,120]; b) increases in intracellular calcium exert a negative feedback on NMDA channels [121]; c) copper modulates the affinity of hippocampal NMDARs for the co-agonist glycine pocampus [122]. As a result, changes in the balance between these trace elements and ions may affect the NMDAR and NMDARmediated neuroplasticity and excitotoxicity, which may be linked to schizophrenia.

Conclusion
All in all, the pathogenic IL-6/IL-23/Th17 axis, which may be induced by microbiota, contributes to neurocognitive deficits and the phenome of schizophrenia, especially MNP, due to its key role in peripheral inflammation, gut and BBB permeability, neuroinflammation and ensuing neurotoxic effects on CNS circuits. Such effects will appear or be more prominent in subjects with CIRS deficits, including lowered IL-10, magnesium, calcium, and zinc.